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Abstract 

It is shown that certain sum rule identities exist which relate correlation func- 
tions for n Potts spins on the boundary of a planar lattice for n > 4. Explicit 
expressions of the identities are obtained for n = 4,5. It is also shown that 
the identities provide the missing link needed for a complete determination of 
the duality relation for the n-point correlation function. The n = 4 duality 
relation is obtained explicitly. More generally we deduce the number of cor- 
relation identities for any n as well as an inversion relation and a conjecture 
on the general form of the duality relation. 
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The Potts model JT], which is a generalization of the two-component Ising model to q 
components for arbitrary q, has been a subject matter of intense interest in many fields 
ranging from condensed matter to high-energy physics. For reviews on the Potts model 
and its relevance see, for example, |2],|3|]. However, exact results on the Potts model have 
proven to be extremely elusive. Rigorous results known to this date are limited, and include 
essentially only a closed-form evaluation of its free energy for q = 2, the Ising model [f§], 
and critical properties for the square, triangular and honeycomb lattices f^-0j. Much less is 
known about its correlation functions. 

In this Letter we report on new sum rule identities for the Potts n-point correlation 
function. Specifically, we show that, as a consequence of being a many-component system, 
the correlation functions of Potts spins on the boundary of a planar lattice must necessarily 
satisfy certain identities when n > 4. We further show that these identities lead to the 
complete determination of a correlation duality relation which, in its simplest form, has 
proven to be useful in determining the equilibrium crystal shape of the Ising model ]8]|| . Our 
results are very general and hold for any planar lattice or graph with arbitrary interactions. 

Consider the g-state Potts model on a planar lattice, or graph, C of N sites and E edges. 
Let i,j, ■ ■ ■ , m, £ be n sites on the boundary ordered as shown in Fig. 1, and let crj denote 
the state of the spin at site i. Two spins of £ at sites i' and j' interact with an interaction 



Kij5(fTi/, ay), where a[, a'j = 1, 2, q. Define the n-point correlation function [T(| 

P n (a, a',---, aW) =< 6{a u a)5{a j , a') ■ ■ ■ S(a e , o») > (1) 

as the probability that the n spins are in respective definite spin states a, a', ■ ■ ■ a^ n \ In 
particular, the correlation function 

T n = q n P n (cr,a,---,a)-l (2) 

vanishes identically if the n spins are completely uncorrelated. 

It is convenient to write P%j-t = P n (h j, • • • £) = Zij-t/Z, where i,j,---,£ = 1, 2, • • • , q, 
Z is the partition function, and Zij...g the partial partition function, namely, the sum of 



Boltzmann factors with the boundary spin states fixed at i,j, ■■-,£. Then we have the 
following theorem. 

Theorem: (i) The correlation functions P n , n > 4, are related by certain sum rule identities. 
Particularly, for n = 4 and 5, the identities are 

-P1212 — -P1213 + -P2131 — -P1234 (3) 
-P12112 = -P12113 + -P21331 — -P12334 

-P12123 — -P12134 + -P21314 — -P12345 (4) 

and eight other relations obtained by cyclically permuting the five indices in (d). 

(ii) The number of correlation identities for a given n is a n = b n — c n , where b n and c n 
are generated respectively from the generating functions 

00 



exp(e* - 1) = E b nt n /n\ (5) 

n=0 

00 

(1 _ v^34t)/2t = Cn t n . (6) 

ra=0 

Proof. The identity (|) is equivalent to 

-^1212 = Z 12 13 + Z 2 131 — ^1234; (7) 

which we represent graphically in Fig. 2. Consider the high-temperature expansion of 
in the form [|ll|] of 

^=E? n(G) n(^'-i)- (8) 

G i',j> V 7 

Here, as a consequence of the fact that the four boundary sites are fixed in definite spin 
states, the summation is taken over all graphs G C £ in which there are n(G) clusters 
excluding those connected to the four boundary sites. 

Apply the expansion (||) to the four Z's in (0). It is clear that, as a consequence of 
C being planar, we have ^1212 — T\ + T-i + T 3 where 7\ is the sum of graphs where sites 
i and k belong to the same cluster, T 2 those graphs where sites j and £ belong to the 
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same cluster, and T 3 graphs i,j, k,£ all belong to different clusters. It is also clear that we 
have Z1213 = T x + T 3 , Z2131 = T 2 + T 3 and Zi 234 = T 3 . The identity (0), and hence (||), 
now follows as a sum rule condition. Clearly, the existence of (|3|) is a consequence of the 
planar connectivity topology. It can also be shown that all n = 5 identities are generated 
by inserting one boundary site to the diagrams in Fig. 2, resulting in identities (^) shown 
graphically in Fig. 3. One can proceed in a similar fashion to derive sum rules for n > 6, 
and thus we have established (i). We remark that the sum rules manifest themselves only 
for q > 4, and therefore do not apply to the Ising model. 

To enumerate a n , the number of correlation identities for a given n, it is instructive to 
consider the case n = 4. First, by enumeration we find that there are 15 distinct Z^i. For 
each ZijM we construct its graph as in Fig. 2 and connect sites in the same state by drawing 
connecting lines exterior to C, resulting in a "connectivity" of the four points. (There is no 
distinction in connectivity topology between drawing connecting lines within or exterior to 
C). A well-nested connectivity, or well- nested Z for brevity, is one in which the connecting 
lines do not intersect [12"]. For n = 4, 14 of the 15 Z's, which are shown in Fig. 4, are 



well-nested. Only Z1212 which, for precisely the same planar topology reason noted in the 
above, is not well-nested. 

More generally for a given n-point correlation function Zij... m g, or Z for brevity, one 
connects in its graph sites in the same state to arrive at an n-point connectivity. Let there 
be altogether b n distinct connectivities of which c n are well-nested. To each Z which is not 
well-nested, we follow the procedure describe in the above, namely, expanding graphically 
in a high-temperature series. Since all graphs in the expansion do not contain intersecting 



lines, by applying the principle of inclusion-exclusion |T3[ we eventually arrive at a sum 
rule expressing the particular correlation function in question in terms of well-nested ones. 
This gives rise to an identity for this particular Z. Furthermore, since each Z has a unique 
graphical expansion, all identities are distinct. It follows that the number of sum rule 
identities, a n , is equal to the number of Z's which are not well-nested, namely, b n — c n . 
The number c n has been evaluated by Blote and Nightingale [|l2j in a consideration of 



the transfer matrix formulation of the Potts model, and is found to be generated by (|6|). To 
enumerate b n we note that it is precisely the number of ways that n objects can be partitioned 
into indistinguishable parts. Let there be m v parts of v objects each, v = 1,2, • ■ -. Then 
we have b n = J2m„=o ' IY^=i[{ n ^)/{ l/ ^) m ' J ' m ^} where the prime over the summation indicates 
the condition Y^™=i vrriv = n. This leads to the generating function (]5[) and thus establishes 
(ii). Particularly, we find a 4 = 15 - 14 = 1, a 5 = 52 - 42 = 10, a 6 = 203 - 132 = 71, a 7 = 
877 — 429 = 448. We have verified these numbers by explicitly enumerating all connectivities 
for n < 6. 

Duality relation for P n : 

It has been known for some time that the two-point correlation function of an Ising 
model is related to its counterpart in the dual space. The usual derivation of this relation 
involves embedding expansions of the correlation functions on the lattice followed by an 
explicit term- by-term identification [14J15]. In a recent paper one of us |I(J introduced a 



new approach to this problem which invokes only a repeated use of an elementary duality 
consideration ||16|| . The new approach, which is very general, also permits the extension 



of the duality analysis to the Potts model for n = 2, 3 ||10|| . However, an extension of the 
analysis of [1(| to n > 4 ran into an apparent snag of inadequacy of conditions [|17]]. Here 
we show that the correlation identities derived in the above provide the missing link, and 
with the help of these identities we determine the duality relation for any n. 



The consideration of [pi is based on the fundamental duality relation |L6 



Z = qCZ* (9) 

relating the partition function Z of any planar lattice, or graph, to the partition function 
Z* on the dual. Here, C = q~ N * fledges ( eJ ^ 3 ~ 1); with iV* being the number of sites of the 
dual and the product taken over all edges. The interaction K*j dual to is given by 

(e*« - l)(e*« - 1) = g . (10) 

Starting from £ we consider a lattice C* formed by introducing n spins a, f3, j ■ ■ ■ , 5 to 



the boundary of the dual of £ (Cf. Fig. 1), each interacting with neighboring dual spins 
within C. (Note that C* now has iV* + n — 1 sites and is not the dual of C.) Let Z*p s be 
the partial dual partition function of C* with the n boundary spins fixed in the respective 
definite states. Our goal is to obtain a duality relation in the form of a linear transformation 
relating the %... m £ to Z* aM ... s . 

Regard the b n well-nested connectivities (such as those shown in Fig. 4 for n = 4) 
as auxiliary lattices, and apply the fundamental duality relation to each one of them |17| . 



Applying the duality on L itself, for example, we obtain (§) which can be written as an 



equation relating linear combinations of the Z and Z* ||10|| . Applying the duality to the 
well-nested connectivity C n in which all n points are connected to a common point with 
interactions K as in £4 shown in Fig. 4(b), we obtain 

^aux(n)= (^)(«-ir^ux(n)> (H) 

where u = e K , and Z aux ( n ) and Z* ux ( n ) are respectively the partition functions of C n and 
its dual. Now, both sides of (|TlD are polynomials of degree n in u. Since ([TT|) holds for 
arbitrary u, the coefficients of all powers of u must be equal. However, it suffices to equate 
only the coefficients of the highest power of u (equating other coefficients leads simply to 
linear combinations of equations to be obtained from other connectivities). On the LHS we 
have ^ a ux(n) = q(u + q — V) n Z\\...\ + 0{u n ~ l \ and find the coefficient of u n to be qZ\\_x- 
On the RHS we have (u - l) n Z* ux(n ) = q(u + q — l) n Z* 1 ... l + other terms, after using fllOD, 
and the coefficient of u n is a linear combination of the b n Z*. This leads immediately to an 
expression for Z X \...\- For n = 4, for example, we obtain 
C 



Znn — -5- 



Z*U1 + Ql(^21U + ^1211 + ^1121 + ^1112) + ?l(^1122 + ^122l) 

+qiq 2 {Z* ll2Z + 

^2113 + ^2311 + ^123l) + 9l(^1212 + QzZ^lZ + ^2^2131 + ?2?3^1234) 

= {1 + qi (l, 1, 1, 1) + q x (l, 1) + qi q 2 (l, 1, 1, 1) + qi (l, q 2 , q 2 , q 2 q 3 )}, (12) 

where q m = q — m,m = 1,2, • • •, and in the last line we have introduced a short-handed 
notation. An immediate consequence of ( |i~2D is the result 

6 



r 4 = gi(p2iu + P1211 + P1121 + P1112 + P1212 + P1122 + P1221) 

+9l?2(Pll23 + P2113 + P2311 + P1231 + P1213 + P213l) + 9l9293£>1234, (13) 

where we have introduced (|S|), = qZ^ ul , as well as p a/ 3 7 s = Z*p S /Z* 1U . For general n 
the consideration of C n leads to 

r n = qi Yl P211-1 P231-1 H 1" ?1?2 • • • 9n-lPl23-n ( 14 ) 

where the meaning of the summations is obvious. 



Applying (1C) to all c n auxiliary lattices of well- nested connectivities in this fashion and 
equating the coefficients of the highest power of u in each case, we obtain c n equations for 
the b n unknown Z's. Since c n < b n for n > 4, it appears that there are more unknowns 
than equations and that the equations are inadequate [HJ. However, after combining the 



c n equations with the b n — c n sum rule identities, we have precisely b n equations, and the 
duality relation can now be completely determined! 

In the case of n = 4, the solution of the 15 equations leads to, in addition to (p"2]), 

2*2111 = {1 + (-1, -1, 9l, ?l) + (<?1, -1) + "I, -1, -1) - (1, 92, 92, 92<?3)} 

^1122 = {1 + (-1, 9i, -1, 9i ) - (1, 1) - 92(1, 1, 1, 1) + (91, qiQz, -92, -9293)} 

^1123 = {1 - (1, -9i, 1, 1) - (1, 1) + (2, 2, -92, -92) + (-1, -92, 2, 2q 3 )} 

%vm = {1 - (1, 1, 1, 1) + 9i(l, 1) - 92(1, 1, 1, 1) + 9i92, 9i92, 9293?")} 

£1213 = {1 - (1, 1, 1, 1) + (-1, 91) + (2, -92, 2, -92) + Q{qi, s, s, q 3 t)} 

Znzi = {1 - (1, 1, 1, 1) + (91, -1) + (-92, 2, -92, 2) + Q{q u s, s, q 3 t)} 

Z l23 4 = {1 - (1, 1, 1, 1) - (1, 1) + 2(1, 1, 1, 1) + Q[r, t, t, 93(2 - 59)]} (15) 

where Q = l/(q 2 — 3q + 1), r = 2q — 1, s = q 2 — 4q + 2, t = q 2 — 59 + 2. Expressions for 
{Z 1211 ,Z 1121 ,Z 1U2 },{Z 1221 } and {^2113,^2311,^1231} are given by cyclic permutations. We 
remark that a closer examination shows that all Z^l except the last four in ( ff3J ) can be 
determined without the use of the identity (|7|). 
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The solution (|P2|) and (|T5|) can be written more compactly by using the fact that Z*, as 
partial partition functions of £*, satisfy the same sum rules as the Z. Particularly, we have 
Z* 212 = Z\ 2Vi + Z 2131 — ^1234- Using this relation and rewriting ([]]) as 

Pa{(JI, 02, CT3, 04) = ^1234 + ^1123^12 + ^2113^23 + ^2311^34 + ^1231^14 
+y4l213<5i3 + ^2131(524 + -4ll22012034 + ^4l221^14^23 + Al212^13^24 

+^2111^234 + ^1211^134 + ^1121^124 + ^1112^123 + ^1111^1234 (16) 

where 5 i2 = S(a%, 02)? ^123 = ^12^23, etc. we obtain after some algebra 
-4i234 = <T 4 [1 - (P2111 + P1211 + P1121 +P1112 + P1122 + P1221 + P1212) 



+2(pil23 + Pl231 + P2113 + P2311 + Pl213 + P213l) - 6^1234] 



^4ll23 


= q~ 


~ 3 (pi211 


- Pl231 - 


P2311 - P1213 + 2^1234) 




^2113 


= q~ 


~ 3 (pil21 


- P1123 - 


Pl231 - P2131 + 2^1234) 




^2311 


= q~ 


~ 3 (Plll2 


- P2113 - 


P1123 - P1213 + 2^1234) 




^1231 


= q~ 


~ 3 G»2111 


- P2311 - 


P2113 - P2131 + 2^1234) 




^1213 


= q~ 


~ 3 (Pl221 


- P1231 - 


P2113 +P1234), -42i3i = q~ 


~ 3 G»1122 - P1123 - P2311 + P1234) 


A\\22 


= q~ 


~ 2 (Pl213 


-P1234), 


Arm = q~ 2 (P2131 -Pl234^ 


), A1212 = 


^2111 


= q~ 


~ 2 (Pll23 


-P1234), 


^1211 = <? _2 (P2113 -Pl234^ 


) 


^4ll21 


= Q~ 


" 2 (P2311 


-P1234), 


^1112 = q~ 2 (Pl231 -Pl234^ 


), A 1U1 = g _1 pi234- (17) 



For general n we write in analogous to flTop 

Pn(&l, &2, ■ ■ ■ , 0"n) = A 12 ... n + Au23...(n-l)^12 H h ^ll-l<W-n (18) 

and similarly for boundary spins a a , erg, ■ ■ • , as of £*. 

P*((Ja, 0/3, • • • , °s) = A* 12 ... n + ^1123-(n-l)<W H ^ A* n ... v 5 a p...5- (19) 

Regard the diagram in Fig. 1 as representing Ay...^ and construct for each A the associated 
connectivity as previously described. Then we are led to the following working conjecture: 

8 



Conjecture: 

Aij...£ = cL l ~ m A* a p,.. s if the connecitivity is well — nested, 

= 0, otherwise, (20) 

where m is the number of distinct indices in {i,j, ■■■,£}. The conjecture is readily verified 
for n = 2,3,4. In practice, for any given n, one can solve A*o 5 from ([190 in terms of 



P<xB-5 = Q Pa/3-s by a PPly m g the principle of inclusion-exclusion jl3| . Details will not be 
given. 

An inversion relation: Since a(3 r j ■ ■ ■ 5 are boundary sites of £*, the transformation relating 
Z* to Z, an inversion process, is given precisely by the same transformation relating Z to 
Z*. Now C* has N* + n — 1 sites and its dual has N — n + 1 sites. Also C and L* have the 
same number of edges. Therefore, we have 



/T\( e K ij — V 

Zm---rne= ry* Hn J ) E M(ij...ml\a(3 1 ...5)Z* aPT .. s , (21) 

V " / af3~f-5 



Z U r -s = ( 3-n^ +in -l ) E M(a/3 7 • • • 6\tij ■ ■ ■m)Z eij ... m , (22) 

where M is a b n x h n matrix. From Fig. 1 we observe that the resulting spin indices on £ after 
the inversion is a counter-clockwise cyclic permutation of the original ordering. Substituting 
(f22|) into (^T|) and making use of ( |T0"D and the Euler relation E = N + iV* — 2, we are led to 
the identity 

M 2 (ijk ■ ■ ■ £\i'j'k' ■■■£>)= f^SifSjv ■ ■ ■ <W, (23) 

which we refer to as an inversion relation. We have explicitly verified this inversion relation 
for n = 2, 3, 4. It can also be shown that the n = 4 inversion relation can be used to deduce 
the sum rule identity (0). 
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Figure captions 

Fig. 1. A planar lattice C and n sites i,j, ... m, £ on the boundary. 
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Fig. 2. Graphical representation of the sum rule identity (||). 
Fig. 3. Graphical representation of the sum rule identities (f|). 

Fig. 4. The 14 well-nested connectivities for n = 4 corresponding to (a) ^1234, (b) Z\m, 
(c) Z U12 occurring 4 times, (d) Z U2 3 occurring 4 times, (e) Z U i3 occurring 2 times, and (f) 
Zn22 occurring 2 times. 
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